Palatini formulation of modified gravity with a nonminimal curvature-matter coupling 
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We derive the field equations and the equations of motion for scalar fields and massive test particles 
in modified theories of gravity with an arbitrary coupling between geometry and matter by using the 
Palatini formalism. We show that the independent connection can be expressed as the Levi-Civita 
connection of an auxiliary, matter Lagrangian dependent metric, which is related with the physical 
metric by means of a conformal transformation. Similarly to the metric case, the field equations 
impose the non-conservation of the energy-momentum tensor. We derive the explicit form of the 
equations of motion for massive test particles in the case of a perfect fluid, and the expression of the 
extra-force is obtained in terms of the matter-geometry coupling functions and of their derivatives. 
Generally, the motion is non-geodesic, and the extra force is orthogonal to the four-velocity. It 
is pointed out here that the force is of a different nature than in the metric formalism. We also 
consider the implications of a nonlinear dependence of the action upon the matter lagrangian. 

PACS numbers: 04.50.-|-h,04.20.Cv, 95.35.-|-d 



I. INTRODUCTION 

A promising way to explain the recent observational 
data [l|, 0] of the accelerated expansion of the Universe 
and of dark matter is to assume that at large scales 
Einstein's general theory of relativity breaks down, and 
a more general action describes the gravitational field 
Ky\. Theoretical models in which the standard Einstein- 
Hilbert action is replaced by an arbitrary function of the 
Ricci scalar R, first proposed in [j, have been exten- 
sively investigated lately. For a review of /(i?) gener- 
alized gravity models and on their physical implications 
see Q . The possibility that the galactic dynamic of mas- 
sive test particles can be understood without the need 
for dark matter was also considered in the framework of 
f{R) gravity models p. 

A generalization of the f{R) gravity theories was pro- 
posed in [3] by including in the theory an explicit cou- 
pling of an arbitrary function of the Ricci scalar R with 
the matter Lagrangian density Lmat [8|- As a result 
of the coupling the motion of the massive particles is 
non-geodesic, and an extra force, orthogonal to the four- 
velocity, arises. The connections with Modified Orbital 
Newtonian Dynamics (MOND) and the Pioneer anomaly 
were also explored, and it was suggested that the matter- 
geometry coupling may be responsible for the observed 
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behavior of the galactic rotation curves. In fact, the non- 
minimal curvature-matter coupling has received much at- 
tention lately. For instance, the model was extended 
to the case of the arbitrary couplings in both geome- 
try and matter in [9|. The implications of the non- 
minimal cou pling on the stellar equilibrium were also in- 
vestigated in [lOJ I , where constraints on the coupling were 
obtained. An inequality which expresses a necessary and 
sufficient condition to avoid the Dolgov-Kawasaki insta- 
bility for the model was derived in [llj. The relation 
between the model with geometry-matter coupling and 
ordinary scalar-tensor gravity, or scalar-tensor theories 
which include non-standard couplings between the scalar 
and matter was studied in [ij] • In the specific case where 
both the action and the coupling are linear in R the ac- 
tion leads to a theory of gravity which includes higher 
order derivatives of the matter fields without introduc- 
ing more dynamics in the gravity sector [l3| . The equiv- 
alence between a scalar theory and the model with the 
non- minimal coupling of the scalar curvature and matter 
was considered in [ij]. This equivalence allows for the 
calculation of the Parameterized Post- Newtonian (PPN) 
parameters /3 and 7, which may lead to a better under- 
standing of the weak-field limit of f{R) theories. 

The equations of motion of test bodies in the nonmin- 
imal coupling model by means of a multipole method 
were also derived in [151]. Furthermore, the energy con- 
ditions and the stability of the model under the Dolgov- 
Kawasaki criterion were studied in [16|. The perturba- 
tion equation of matter on subhorizon scales in models 
with an arbitrary matter-geometry coupling was used to 



constrain the theory from growth factor and weak lensing 
observations [13] • In particular, the age of the oldest star 
clusters and the primordial nucleosynthesis bounds were 
used in order to constrain the parameters of a toy model. 
The possibility that the behavior of the rotational veloc- 
ities of test particles gravitating around galaxies can be 
explained in the framework of modified gravity models 
with non-minimal matter-geometry coupling was consid- 
ered in [1^ . Analogous nonlinear gravitational couplings 
with a matter Lagrangian were also considered in the 
context of proposals to address the cosmic accelerated 
expansion [19|. 

The subtle issue of a correct definition for the mat- 
ter Lagrangian of the theory and of the definition of the 
energy-momentum tensor in the presence of a nonmini- 
mal curvature-matter coupling needs to be pointed out. 
In a recent paper [20|, the authors argued that a "natu- 
ral choice" for the matter Lagrangian density for perfect 
fiuids is Lmat — P, based on Refs. [2l|,|23], where p is the 
pressure. This choice has a particularly interesting ap- 
plication in the analysis of the curvature- matter coupling 
for perfect fluids, which imposes the vanishing of the ex- 
tra force p, [23] . Despite the fact that Lmat — P does 
indeed reproduce the perfect fluid equation of state, it is 
not unique |24| . Another choice includes, for instance, 
^m = —p 123 . [25I |26| , where p is the energy density (see 
Ref. [22, [2J] for details). For a review of modified f{R) 
gravity with geometry-matter coupling see [231 . 

In the literature different approaches in f{R) modified 
theories of gravity are used. These include the metric 
formalism, where the action is varied with respect to the 
metric; the Palatini formalism, where the metric and the 
connections are treated as separate variables, and the 
Lagrangian is varied with respect to both to derive the 
field equations; and the metric-affine formalism, which 
generalizes the Palatini variation, where the matter part 
of the action depends and is varied with respect to the 
connection. In this work, we use the Palatini formal- 
ism, that treats the metric and the affine connection as 
independent geometrical quantities. When the Einstein- 
Hilbert action is used, the Palatini variational principle 
leads to the Einstein equations, as in the standard metric 
variation. This is not true, however, for a more general 
action. When used together with an f{R) Lagrangian, 
the Palatini formalism leads to second order differential 
equations instead of the fourth order ones that one gets 
with the metric variation. At the same time, in vacuum, 
they straightforwardly reduce to standard General Rela- 
tivity (GR) plus a cosmological constant. This ensures us 
that, firstly, the theory passes the solar system tests, and 
secondly, that interesting aspects of GR like static black 
holes and gravitational waves are still present. Thus 
there is no real criterion so far about which one of them is 
better to use. Additionally, the Palatini variation seems 
to be more general since it yields GR without the need to 
specify the relation between the metric and the connec- 
tion. The Palatini formalism for f{R) modified gravity 
models has been intensively investigated recently [28] . 



As alternatives to dark energy, these models run into 
problems with structure formation [2^. While the mod- 
ified gravity sector can support an accelerating universe 
without dark energy, the cosmological perturbations will 
then behave differently than in GR with a cosmological 
constant. In order to comply with the observations of 
the large-scale structure, one is then forced to to fine- 
tune the gravity lagrangian to be essentially GR with 
a cosmological constant. However, it has been shown 
that this can be avoided in the presence of generalized 
dark matter [30|. In particular, the evolution of cosmo- 
logical perturbations can be viable if dark matter had 
effective pressures related to the modified gravity sector. 
Such effect ensues from a nonminimal coupling to curva- 
ture, motivating to study such a coupling particularly in 
the Palatini formalism. In fact, dark matter itself could 
emerge from a nonminimal coupling of baryons to curva- 
ture at supergalactic scales. 

It is the purpose of the present paper to derive the 
gravitational field equations of the generalized f{R) type 
gravity models with non-minimal coupling between mat- 
ter and geometry, which depends on two arbitrary func- 
tions of the Ricci scalar and of the matter Lagrangian, 
respectively. By taking two independent variations with 
respect to the metric and the connection separately of the 
gravitational action, we obtain the field equations and the 
connection associated to the Ricci tensor, which, due to 
the matter-geometry coupling, is also a function of the 
matter Lagrangian. The metric that defines our indepen- 
dent connection is conformally related to the spacetime 
metric, with the conformal factor a function of the matter 
Lagrangian and of the Ricci scalar. Once the conformal 
factor is known, the field equations can be obtained eas- 
ily in both metrics. By taking the divergence of the field 
equations it follows that the energy-momentum tensor of 
the matter is not conserved, and, similarly to the metric 
case, due to the matter-geometry coupling an extra force 
arises. 

The present paper is organized as follows. In Section 
im the field equations of the modified gravity model with 
arbitrary matter-geometry coupling are obtained by us- 
ing the Palatini formalism. The equation of motion of 
test bodies is derived in Section lllll We discuss and con- 
clude our results in Section [IVI In the present paper we 
use the natural system of units with SttG = c = 1 . 



II. FIELD EQUATIONS IN PALATINI 

MODIFIED GRAVITY WITH AN ARBITRARY 

COUPLING BETWEEN MATTER AND 

GEOMETRY 

By assuming an arbitrary coupling between matter and 
geometry, the action for the modified theory of gravity 
takes the form [I] 
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where fi{R) (with i = 1, 2) are arbitrary functions of the 
Ricci scalar R — g'^'^R^i,, while G {Lmat) is an arbitrary 
function of the matter Lagrangian density Lmat- The 
matter Lagrangian is a function of the metric tensor g 
and of the physical fields -07 while the Ricci tensor R^^ 
is expressed solely in terms of the connection T. The 
only requirement for the functions fi{R), i = 1,2, and 
G, is to be analytical functions of the Ricci scalar R and 
of the matter Lagrangian density Lmat , respectively, that 
is, they must possess a Taylor series expansion about any 
point @, |31[ . The Ricci tensor is defined as p, |32| 
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with the connection F^^, obtained through the indepen- 
dent variation of the gravitational field action given by 
Eq. (H)), and not directly constructed from the met- 
ric by using the Levi-Civita prescription. The energy- 
momentum tensor of the matter is introduced as 



T — — - 
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thus giving 
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The Palatini formalism consists in taking separately 
two independent variations with respect to the metric 
and the connection, respectively. The action is formally 
the same, but the Riemann tensor and the Ricci tensor 
are constructed with the independent connection. By 
varying the action ([T} with respect to the metric g we 
obtain 



[![{R)+2f!,{R)G{Lmat)]Rt.. 

\h{R) + h{R)G{Lmat) 
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where the prime denotes a derivative with respect to 
the argument, i.e., fl{R) = dfi{R)/dR and G" {Lmat) — 
dG [Lmat] I dLmat- In the limit /2 (i?) = 1, G {Lmat) = 
Lmat we obtain the Palatini formalism field equations of 
the standard f{R) modified theory of gravity [a, Q- In 
the case fi{R) = R, f2{R) = 1, and G {Lmat) = Lmat, we 
obtain the field equations of standard general relativity. 
The next step in the Palatini formalism requires the 
variation of the action with respect to the connection F. 
The variation can be done by using the identity 



SR^. = Va (<5f^.) - V^ (^f^;,) , 



(6) 



where Va is the covariant derivative associated with the 
connection F. 



By taking the variation of the action ([T]) with respect 
to the connection F we obtain 
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where we have denoted 
'1 



^f,^ 



f[{R) + f^{R)G{Lr^ 



(7) 
(8) 



By integrating by parts we obtain 
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The first term in 5S/5T is a total derivative, and thus 
it can be discarded. Hence the variation of the action 
with respect to the connection gives 



^^{A^^''5l-A^''6li)] =0 



(10) 



Equation (llOp can be further simplified if we take into 
account that when a = \ the equation is identically zero. 
Taking A 7^ a, we obtain 

V4V^[/{(i?)+2/^(i?)G(i„,a,)].9^"} =0. (11) 

Equation (|TT|) shows that the connection is compatible 
with a conformal metric. By defining, according to |5|,|8|, 
a new metric h^^, conformal to g^jj, given by 



h^, = [f[{R) + 2f'^{R)G {Lmat)] 5^., (12) 

we obtain 

yfHU^'' = ^—g[f[{R) + 2f^{R)G{Lmat)]g^\ (13) 

where h is the determinant of the metric h^v Thus 
Eq. pT|) becomes the definition of the Levi-Civita con- 
nection F of h^i, , giving 



f ' - ^h^' {d.h 
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With the use of the explicit form of /i^^ we obtain 
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where we have denoted 

F^F{R, Lmat) = f[{R) + 2f^{R)G {Lmat) ■ (16) 
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In terms of the Levi-Civita connection F^^, associated 



r^^ = ^9^'' {di^gpp + dpg^p - dpgp^) , 
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r^^, can be expressed as 

f ^. = r^, + d. In ^5l + d^ In VFSt - g^.g^'d^ In Vf. 

(18) 
The Ricci tensor i?^i, is given in terms of the tensor 
Rfj,i, constructed from the metric by using the Levi-Civita 
connection Eq. ([T7l) by [5| 
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while the Ricci scalar and the Einstein tensor can be 
obtained as 

i? = i?(.9)+3inF+^ij(V^i^)(V'^F), (20) 

and 



G^. = ^p.-^3M-i?-G^.(ff) + ^;^(V^^)(V,i^) 
^\ (V^V. + .g^.D) F - \^g^.. (VaF) (V^F) ,(21) 

respectively, with all the covariant derivatives taken with 
respect to the metric g^J,v. 

With the use of the expression of the Einstein tensor 
given by Eq. (PT|) the gravitational field equation Eq. ([5]) 
can be written as 
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where for notational simplicity, we introduced the nota- 
tion 
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By substituting the expression of the Einstein tensor 
given by Eq. (PT|) into the field equation Eq. (|22p we ob- 
tain the gravitational field equation of modified gravity 
with a nonminimal coupling between matter and geome- 
try in the Palatini formalism as 

2K 11 R\ , f2{R)G{L„,at) ^ 

— 2F°^-^) ^'^ + F ^-- (2^^ 

Note that the field equation Eq. ([5]) can be written as 

FR,,, - ^K {R, Lmat) .9p. = f2iR)G' {Lmat) T^.- (25) 
By taking the trace of Eq. ([^5)) we obtain 

FR - 2K {R, Lmat) = f2{R)G' {L^at) T, (26) 



where T = T^ is the trace of the energy-momentum ten- 
sor. With the use of Eq. ([T^ we find the equation deter- 
mining _R as a function of T and the matter Lagrangian 



^(5) 



sIdf- 



3 J_ 

2F2 



(V^F) (V^F) 



f,! 7^/n r \ f2{R)G' (Lmat) rr, 
-2 — K {R,Lraat) = -^ J ■ 



(27) 



III. EQUATION OF MOTION AND THE 
EXTRA FORCE 



A. Non-geodesic motion 

By taking the covariant divergence of Eq. (I^5t we ob- 
tain 

(V^F) R^, + FV^'R^, ~ \g^.uV^K {R, Lmat) = 

= [V^/2(i?)G" {Lraat)] T^. + f2{R)G' (L^at) VT^,. 

(28) 
With the use of the identity 8] 
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(29) 



we can write 



FWR^, + (V^F) R^, = -g^^FV^R. (30) 



Using this in Eg. ([^5)) we then obtain that 
^A^n = TTm^r -A~{^^.h{R)G'{Lraat))T'^. 

+ L,natG'{Lraat)^.h{R)-h{R)V,gY (31) 

where for notational convenience we introduced the 
shorthand Q = G{Lmat) - L^atG' {Lmat), which vanishes 
when the action is linear in the matter lagrangian. 

For a perfect fiuid, described by an energy density 
e and a thermodynamic pressure p, the matter energy- 
momentum tensor is given by 



rlr^ = (e + P) "m"'' + Wm'" 



(32) 



where the four- velocity, u^, satisfies the conditions 
u^u^ = — 1 and u^^u^-i, = 0. The trace of the energy- 
momentum tensor is given by T = 3p — e. We also in- 
troduce the projection operator /i^a = 5/^a + u^ua from 
which one obtains h^\u^^ — and T^^'^h^x = ph\^ respec- 
tively. 

By contracting Eq. ([^5]) with the projection operator 
hf^x, one deduces the following expression 

[e + p] g^xu'V.u^ ~ i^uPWx - u'ux) = 

^ (5^ - U-'ux) X {pV, (/2(i?)G' {L,nat)) 
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(33) 



Finally, contraction with g"^ gives rise to the equation 
of motion for a fluid element 



ds 



du" 
ds 
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(34) 



where we have introduced the space-time connection F" , 
which is expressed in terms of the Christoffel symbols 
constructed from the metric g^„, and where the extra- 
force is defined as (dropping the explicit dependence of 
G upon Lmat)- 



r 



h" 



e + p 
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As one can see from Eq. ([55]) . the extra force /" is or- 
thogonal to the four- velocity of the particle, 



ru^ = 



(36) 



The extra force-four velocity orthogonality relation fol- 
lows from the properties of the projection operator. This 
is consistent with the usual interpretation of the force, 
according to which only the component of the four-force 
that is orthogonal to the particle's four-velocity can in- 
fluence its trajectory. The presence of the extra force /" 
implies that the motion of the particle is non-geodesic. 
For /" = we recover the geodesic equation of motion. 
The usual gravitational effects, due to the presence of an 
arbitrary mass distribution, are assumed to be contained 
in the term a^ = F^ow"u^. 

We can now separate three contributions to the force. 
The first term in ([35]) is the usual which is due to the 
gradient of pressure. The second term is due to the 
nonminimal curvature coupling to gravity and of course 
vanishes in general relativity. We note that the result 
reached in the metric formalism applies also here, that 
when L^at = Pi the extra force in fact disappears. Inter- 
estingly, this is the case in the usual lagrangian descrip- 
tion of scalar fields and can hold for more general perfect 
fluids as well, depending upon the chosen prescription. 

Furthermore, the two following terms, consisting of the 
last line in ([55]) . are due to the new self- interactions of 
matter and disappear when the action is linear in the 
matter lagrangian L^at- AH of the three terms can be 
considered as contact forces. For a test particle, these 
vanish when in vacuum. This reflects the fact that the 
Palatini theory of gravity can be understood as a mod- 
ified response to matter sources whereas it reduces to 
general relativity (with a cosmological constant) in vac- 
uum. In other words, though the geodesies of particles 
are modified within a matter distribution, there are no 
new propagating degrees of freedom which would medi- 
ate interactions in vacuum and the extra force is absent 
there. 



B. Scalar field 

Let us as an example consider a general scalar field 
theory given by a lagrangian depending upon the field (f). 
The Euler-Lagrange equations of motion read then 
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dLr. 
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The left hand side represents now the extra force. This 
is due to nonminimal curvature coupling if f2{R) is 
not a constant, and due to matter self interactions, if 
G'{L,riat) is not a constant. Consider the purely latter 
case. As found above, now the extra-force entering the 
geodesic equation vanishes, but however there are addi- 
tional terms in the Klein-Gordon equation for the field 
in curved space. Thus the additional terms do not affect 
the motion of a scalar particle. On the other hand, in 
the latter case, we obviously obtain new effects already 
in Minkowski space, and these can be regarded as non- 
gravitational extra force. 

As an explicit example, we let us consider an ac- 
tion given by an exponential of the Einstein-Hilbert la- 
grangian. This we fix the exponential dependence of the 
nonminimally coupled matter sector as 
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where A is a suitable energy scale. Let us keep the form 
of the scalar field lagrangian as in the canonical theory. 



Lraat = \ {d^f - V {c^) . 

The Klein-Gordon equation then becomes 



(40) 



T^d.Rd^cl) 



= 0. (41) 



The corrections appear within the square brackets and 
are suppressed by the scale A. The propagator is mod- 
ified in a nonlinear way. Note that only the last term 
in the second line drops when we consider fiat instead of 
curved space. 



IV. DISCUSSIONS AND FINAL REMARKS 

In the present paper we have considered the gravita- 
tional field equations for a generalized f{R) type grav- 
ity model with a geometry-matter coupling in the frame- 
work of the Palatini formalism. Similarly to the metric 



case [3, Q the energy-momentum tensor of the matter 
is not conserved. Palatini type f{R) theories have spe- 
cial properties that make them especially interesting for 
addressing strong gravity phenomena such as the early 
Universe or stellar collapse processes [2^. The coupling 
of the matter Lagrangian with the curvature term leads 
to an extra- force term in the equation of motion of mas- 
sive particles in gravitational fields, which essentially de- 
pends on the allowed functional forms for the geometry- 
matter coupling. If, for example, the geometry-matter 
coupling could generate some repulsive forces, then one 
could construct non-singular cosmological models that 
are non-singular at high curvatures, or even non-singular 
collapsing stars. No new degrees of freedom in the grav- 
itational side or in the matter side (exotic sources) are 
needed to get such repulsive gravitational forces, since 
the extra-force is induced by the coupling between mat- 
ter and geometry. Our results suggest that Palatini type 
theories might play a relevant role in the phenomenology 
of gravitation at both high energies (densities), as well 
as in the very low density limit. On the other hand, the 
assumption of independence between metric and connec- 
tion in the variational process is essential to get second 
order equations for the metric. It is thus reasonable to ex- 
pect that effective descriptions of the gravitational force 
at large/small scales could come in the form of Palatini 
theories. 

In three dimensions and in the Newtonian limit, 
Eq. (j35|) can be formally represented as a three-vector 
equation of the form a = a^ + f where a is the total ac- 
celeration of the particle, Sn is the gravitational acceler- 
ation and / is the acceleration (per unit mass) due to the 
jresence of the extra force. Similarly to the metric case 
3, a Modified Newtonian Orbital Dynamics (MOND) 
33| type acceleration is induced by the presence of the 
geometry-matter coupling. This shows that one possi- 
bility of testing the effects of the matter-geometry cou- 
pling could be in the region of very small accelerations 
of the order of 10~^" m/s^. A supplementary MOND- 



type acceleration can explain the observed behavior of 
test particles gravitating around galaxies. 

As a possible physical application of the Palatini for- 
malism gravitational field equations one could consider 
an alternative view to the dark matter problem, namely, 
the possibility that the galactic rotation curves, and the 
mass discrepancy in galaxies and clusters of galaxies, 
which is usually explained by postulating the existence 
of dark matter, can be explained in gravitational models 
in which there is a non-minimal coupling between mat- 
ter and geometry J3, [y, UM ■ Similarly to the metric case, 
the extra-terms in the gravitational field equations of the 
Palatini formalism modify the equations of motion of test 
particles, and induce a supplementary gravitational inter- 
action, which can account for the observed behavior of 
the galactic rotation curves. Therefore, the analysis of 
the motion of the test particles at a galactic scale could 
open the possibility of directly testing the modified grav- 
ity models with non-minimal coupling between matter 
and geometry in the Palatini formalism by using direct 
astronomical and astrophysical observations at the galac- 
tic or extra-galactic scale. In this paper we have provided 
some basic theoretical tools necessary for the in depth 
comparison of the predictions of this model and of the 
observational results. 
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